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Abstract 

In this paper, a new triangular decomposition algorithm is proposed for ordinary differ¬ 
ential polynomial systems, which has triple exponential computational complexity. The 
key idea is to eliminate one algebraic variable from a set of polynomials in one step using 
the theory of multivariate resultant. This seems to be the first differential triangular 
decomposition algorithm with elementary computation complexity. 
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1 Introduction 

A basic problem in computer algebra is to properly represent the solutions for a set of 
algebraic or differential equations and the triangular set is one of the basic ways to do that. 
Let f i ..., f s be polynomials in variables xi,... ,x n . Then it is possible to compute triangular 
sets 7i,...,T r such that 


Zero(/i ..., / s ) = U[ = ! Zero (sat (7p) 

where sat(7i) is the saturation ideal to be defined in section 2 of this paper. Since each % is 
in triangular form, many properties of its solution set can be easily deduced. Triangular de¬ 
compositions also lead to many important applications such as automated theorem proving, 
kinematic analysis of robotics, computer vision, stability analysis of molecular systems, etc. 

The concept of triangular set was introduced by Ritt pjjj in the 1950s and was revived in 
the 1980s by Wu [20] in his work of automated geometry theorem proving. One of the major 
advantage of the triangular decomposition method is that it can be used to give complete 
methods for the radical ideal membership problem of differential and difference polynomial 
ideals, while the well known Grobner basis method does not suit for this purpose. By now, 
various kinds of triangular decomposition algorithms have been proposed for polynomial 
systems nuni ttsiEn], differential polynomial systems [21 El El [T3l [19] , difference polynomial 
systems [9], polynomial systems over finite fields [SOUS], and semi-algebraic sets g]. 

The computational complexity analysis for triangular decomposition algorithms is quite 
difficult and only very limited results are known. For polynomial systems, Gallo and Mishra 
gave a single exponential algorithm to compute the characteristic set for a finitely gener¬ 
ated ideal [6] and Szanto gave a randomized single exponential algorithm to compute the 


1 


triangular decomposition m- The complexity analysis of the commonly used triangular 
decomposition algorithms is not given yet. However, it is shown that if solutions in Z 2 
are considered, then the commonly used triangular decomposition algorithm can be made 
single exponential and practically very efficient [8|. For differential polynomial systems, 
it is generally believed that the commonly used triangular decomposition algorithms have 
non-elementary computational complexity [IHl- 

In this paper, new triangular decomposition algorithms are proposed for polynomial 
and differential polynomial systems. The key idea is to eliminate one algebraic variable 
from a set of polynomials in one step using the theory of multivariate resultant. This 
method was introduced by Yu Grigor’ev to give a quantifier elimination algorithm with 
nice computational complexity HD • In this paper, by adapting this elimination method, we 
give triangular decomposition algorithms for polynomial and ordinary differential polynomial 
systems. In the case of polynomial systems, the algorithm gives an unmixed decomposition 
and has double exponential complexity. In the case of differential polynomial systems, the 
algorithm gives an unmixed radical decomposition and has triple exponential complexity. 
This seems to be the first differential triangular decomposition algorithm with elementary 
computation complexity. 

The rest of this paper is organized as follows. In Section 2, we give a new triangular 
decomposition algorithm for polynomial systems. In Section 3, we give a new triangular de¬ 
composition algorithm for ordinary differential polynomial systems. In Section 4, a summary 
is given. 

2 Decomposition of algebraic polynomial system 

In this section, we give an algorithm which for given polynomials h\, ..., hk € IK[xi,..., x n ], 
gives the decomposition Zero(/ii,..., h &) = U g Zero(sat(^l ( j)), where A q is a regular triangular 
set for each q. Furthermore, the computational complexity of the algorithm is given. 

2.1 Basic definition and property 

Let IK be a field of characteristic 0, and x\ < X 2 < ■ ■ ■ < x n ordered variables. For every 
i € {1,... , n}, we define Kj = IK[xi,..., Xj] to be the ring of multivariate polynomials in the 
variables xi,..., Xj with coefficients in IK. We write deg(/, Xi) for the degree of / in x,;, and 
deg Xi xi 2 ,...,x it (/) f° r the degree of / as the multivariate polynomial in xq,Xj 2 , ... ,Xi t . We 
call the leading variable of /, denoted by lv(/), the greatest variable v G {xi,... ,x n } such 
that deg(/, v) > 0. 

Assuming lv(/) = x*, we call i the class of /, denoted by cls(/). Regarding / as a 
univariate polynomial in IKj_i[xj], we can write / = cxf + r. We call d = deg(/, Xj) the 
leading degree of /, denoted by ldeg(/), and c the initial of /, denoted by ini(/) or If. 

For Let P be a polynomial set and D a polynomial in lK n . For an algebraic closed 
extension held E of IK, let 

Zero(P/L>) = {q G E n | VP G P, P(rj) = 0 A D(rj) / 0}. 

A subset T of IK n is called a triangular set if no element of T lies in IK and for P,Q G T 
with P/Qwe have lv(P) / lv(Q). 

Let T = {Ti,... ,T r } be a triangular set. We always assume lv(Ti) < lv(T 2 ) < • • • < 
lv(T r ). We can rename the variables as iti,... , u q , yi, ..., y r such that q + r = n and lv(Tj) = 
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Hi. ThenT has the following form: 


TiOi,... ,u q ,yi) 

q- _ ( T 2 (u\, •••) Hit 2 / 2 ) ^ 

„ T r (u \,..., Uq, yi ,..., y r ) t 

We call u = {ui ,..., u q } the parameter set of T, and write Ij- = It± ■ ■ ■ It t ■ For a triangular 
set T, the saturation ideal of T is defined to be 

sat(T) = {/ G I„ | 3d € N + , si. if/ € (T)} 

where (T) is the ideal generated by T. 

A triangular set T = [T±,... ,T r \ of form (1) is called regular, if for each 1 < i < r, 
(Ti,..., Tj_i, ini(Tj)) P) K[u] 7 ^ {0} where (T\,..., Tj_i, ini (7))) is the ideal generated by 
Ti,..., Tj_i, ini(Tj) and u is the parameter set of T. 

Lemma 2.1 [7/ For a triangular set T, we have: 

t 

\J sat(T) = P) sat(7j) 

i=l 

where Ti are regular triangular sets having the same parameter set as T, and sat (71) is a 
prime ideal. That is, sat (T) is an unmixed ideal. 

Lemma 2.2 Let T = {Ti,...,T r } 6 e a regular triangular set, u its parameter set, yi the 
leading variable of T, P a polynomial in K[tz, y\,... ,y r \. Then ( P , T) f) IK[«] T {0} if P is 
not identically zero on all irreducible components 0 /Zero (sat (T)) 

t 

Proof: According to Lemma 12.11 ^/sat(T) = f) sat(^). Since P is not identically zero on 

2=1 

all irreducible component of Zero(sat(T)), we have P ^ sat (Qi) for each i. Since sat(^) is 
prime, so ( P , Qf) H K(u) / {0} for each i. Suppose that Q t = (G^i,... , Gi )T ) , then we have 
the following equalities: 


Si,iGi,i + • • • + Si )V Gi t r — A\P + hi(u) 
5 2 ,iG 2 ,i + • • • + S 2 , r G 2 ,r = A 2 p + h 2 (u) 


St,iGt,i + • • • + St )T Gt,r — A t P + ht( u). 

Multiply all the equalities. Since the left hand side of the i-th equality belongs to sat(£?j), 
the product of them belongs to ysat(T). The product of the right hand side is of the form 
AP + h for h = h\h 2 - ■ ■ ht- Then we have h G ( P , ^/sat(T)) P|^[ u ] 7 ^ {0}- Therefore, 
there exists an integer d\ such that ( h) dl G (/, sat(T)). There exists an integer d 2 s.t. 
(ini(Ti)... ini(T r )) d 2 (h) dl G ( P,T ). Since T is regular, we have (ini(Tj), T) f) IK[u] T {0}, 
and then the following equality 

(ini(Ti)... ini(T r ))*F 0 = g( u) + F l T 1 + ... + F r T r 

where g / 0 and g G K[u]. Hence, ( h) dl g G ( P,T ) and ( P,T) P|^[ u ] 7 ^ {0}- □ 

As a consequence of Lemma 12.21 we have: 
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Corollary 2.3 A triangular set T = {Ti,... ,T r } is regular if for each 2 < i < r, ini(T)) is 
not identically zero on all irreducible components of sat (T \,... , T)_i). 

Lemma 2.4 Let T = {T\,... ,T r } be a regular triangular set with parameter set u. If 
M E K[«] . then Zero iff/I t) — ZeroifT / M I-p), where S is the Zariski closure ofS. 

Proof: We first prove the lemma when sat(T) is prime. Introduce a new variable z and let 
I = I-p- We have 

Zero (T/MI) = Zero(T, MIz - 1) D E n . 

So for any / E I(Zero (T/MI)) (1(C) is the ideal of polynomials which vanish on V), f E 
a/ (T, MIz — 1). Let z = jjj, then there exists an integer d such that ( MIf) d E (T), 
so we have ( Mf) d E sat(T). Since sat(T) is prime and M E K[u] so not in sat(T), we 
have / E sat(T). So we have Zero (T/MI) D Zero(sat(T)) = Zero (T/I). It is obvious 
Zero(Ti,... ,T r /MI) C Zero (T/I), so we have Zero (T/MI) = Zero {T/I). Now assuming 

t 

sat (T) is not prime. According to Theorem 1.3 in we have Zero {T/I) = (J Zero {Gi/If), 

2—1 

where Qi is a regular triangular set having the same parameter set with T and sat(^) is 
prime for each i. Then we have 

_ t _ t _ _ 

Zero (T/MI) = (J Zero (ft/I.M) = [j Zero (ft/i*) = Zero (T/I) 

i— 1 2— 1 

and the lemma is proved. □ 

Lemma 2.5 Let T = {Ti,... ,T r } be a regular triangular set. If M is not identically zero 
on all irreducible components of sat(T), then Zero(sat(T)) = Zero(T//r) = Zero(T/M/ 7 -), 
where § is f/ie Zariski closure o/§. 

Proof: It is well known that Zero (sat (T)) = Zero (T/It)- We have Zero (T/Tr) ^ Zero(T/M I-p), 
since Zero(T//r) A Zero(T/MIf). Let u be the parameter set of T, since M is not 
identically zero on all irreducible components of T, according to Lemma 12.21 we have 
(M, T)P|^[ U ] ~f~ {0}- Suppose 

AM + A\T\ + ... + A r T r = H( u) (2) 

Let ^ be a zero of Zero (T/HI-p). Substitute £ into ©, we have M{tf) ^ 0, so ^ is also a 
zero of Zero(T /I-p) and we have Zero(T/ MI-p) D Zero(T/^r)- Since H E K[u], so according 
to Lemma 12.41 we have Zero (T/HIp-) = Zero (T/It)- Therefore, we have Zero(T /M I-p) D 
Zero (T/HI-p) = Zero (T/It)- This completes the proof. □ 

2.2 A quasi GCD algorithm 

We need to use lemma 1 of Em, which is modified slightly to the following form. 

Lemma 2.6 fllf There is an algorithm which for given polynomials hi = ^ffh ijY 3 E K n [T], 
deg 3 , 1 x y (/ij) < d,i = 0,1,... , k, yields such two families of polynomials g q: t E K n , E 
K n [Y] for 1 < (/ < Ai, 0 < t < IV 2 such that 

m 

Zero(/ii,..., h^/ho) = |J Zero(^ 9 , g q ,i ,..., g q ,N 2 /g q ,o) 

q =1 

(J Zero({hjj, i = 0,..., k,j = 0,..., d - l}/h 0 ). 
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Furthermore, we have the following properties. 

(1) . deg( v I' g . Y) > 0, ini('I' ? ) | g qfi . 

( 2 ) - de Sx 1 ,...,x n ,Y(^q)^ e Sx 1 ,...,xn(9q,t) < V{d)]N 1 ,N 2 < kV{d n ) whereV(k) means a poly¬ 
nomial in k. 

(3) . The running time of the algorithm can he bounded by a polynomial in k and d n . 

Now we describe the main steps of this algorithm without proof. One can refer to El 

for more details. Without loss of generality, we assume that degy(hj) > 0 for 0 < i < k. 

d— l _ j 

Since deg a , 1 Xn y{hi) < d, we have hi = hijYT Let hij = h^pY^ and 

3=0 ’ >3=0 

Ui,j = Zero {hi d—i, ■ ■ ■, h±g, —15 • • • ? h‘ 2 ,o, • • • > i; ■ ■ ■ > hij+i/hij') 

for 1 < i < fc, 0 < j < d — 1. Let H = {h\,..., h+ 

Hij {h/.j. ..., h/j} (3) 

Hk+i = {hij, V 1 < i < k and j}. 

Then we have Zero(iL/ho) = Zero(iLfc + i/ho) (J( (J Zero(iLjj//to) fl C/jj). 

l<i<fc,0<j<d-l 

Now we turn to the system Hjj and introduce new variables Iq, Yi to make polynomials 
in Q homogeneous in Y. Yq, Y\. Let 

_ y 

hi = Y$hij(x i,...,x n ,—) 

*0 

hi = YQ deg( ' hl) h l (x 1 ,...,x n ,^ r ),i + l <1 <k 

ho = Y^ eg(ho)+ \^-h 0 (x 1 ,...,x n ,^)-l) 

Y 0 Y 0 

Then ho, hi,... ,hk are homogeneous polynomials in Y, Yq,Y\. The solutions of the following 
homogenous system correspond bijectively to that of Q except ( 1 : 0 : 0 ). 

hi = h i+ 1 = ... = h k = h 0 . (4) 

Here ho, hi,... ,h k are considered as polynomials in Y, Yq,Y\. 

Introduce new variables Uq, U, U± and let h k+ 1 = YqUq + YU + Y\ U\. We rearrange the 
polynomials ho, hi,...,h k w.r.t the degree in Y, Yq, Y\ as go,..., g k -i +2 and 70 > 71 > ... > 
7 fc —«+2 where degyy o yi (. 9 . s ) = j s for 0 < s < k - i + 2 . Since degyy 0 Yl (h k+1 ) = 1 , we can 
assume that g k -i +2 is h k+ 1 . Let 

D = { Y (71 — 1)) + 7o- 

l<Z<min{2,A:—i+1} 

We construct the Macaulay matrix A as the representation of the linear map 


A : TLq © ... © TL k -i+2 hi 


where TLi (respectively TL) is the linear space of homogenous polynomials in Y, Yq, Y\ of degree 
D — 7 ; (respectively D) for 0 < l < k — i + 2, and 

A{fo, ■ ■ ■ , fk-i+ 2 ) = fodO +-1- fk-i+29k-i+2 
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k—i+2 

The matrix A is of size C‘f }+ , 2 x ^ C'| ) _ +2 an d can be represented in a form A = 

1=0 1 

^(num)^y[(for)^ w h ere the elements of the submatrix do not contain U,Uq,Ui. Ac¬ 

tually, A^ num ' ) is the submatrix of A which corresponds to the basis of TLq, ..., TLk~i +1 while 
Aif or ) corresponds to the basis of TLk-i+ 2 - 

About the polynomial system (J4|) , we have the following lemma: 

Lemma 2.7 ,171]/ The rank of the matrix A of the polynomial system 0 is r = Cp +2 - Let 
A be a nonsingular r x r submatrix of A containing rank( A ^ num ^) columns in Afnum)' Then 


D i 

det(A) = Li, where Li = fifiUo + fiU + and c is a constant 

1=1 

where (£j ; o : £* : £i,i) a solution of 0 and the number of occurrences of iifiUo+fiU+fi,iU\ 
in the product coincides with the multiplicity of the solution (£^0 : fi ■ 1 ) of 0- 

To hnd the A in Lemma l2.7l we use a variant of Gaussian algorithm which will compute 
a series of 

W s = {x € : P 1 = ... = P s _i = 0 ,P S ± 0} (5) 

where Pi,P, are polynomials in x, U, Uo, U\ and linearly independent. For x E W s n Uij, 
the determinant 

d 2 

A s = Y J E i i ) U^- i (6) 

i=0 

is what we want. For more details about the variant Gaussian algorithm, one can refer to 

Ha- 

Now we introduce the following quasiprojective varieties 

W® = {x <E W s : E^ = ... = EV-V = 0, E® + 0} (7) 

where Es°\ • • ■, E's ^ are polynomials in x, U, U±. In Ha, it is proved that if we substitute 
Ui = 0, U = —1,Uq = Y into -%r, and denote the polynomial as \I/ S , then for each point 

Es 

x G Ws FI Ui j, the solution of T s as a polynomial in Y is the solution of the polynomial 
system ©• Since the quasiprojective varieties Wp (~l Uij can be divided into a series of 
polynomial systems V t = Zero .. ■, gt,N 2 /9t,o), we have 

Zero (H itj ) D U itj = |J Zero(T t , g t> 1 ,..., gt,N 2 /9t,o)■ 

t 

If dq = 1, we can delete that component and finally obtain the decomposition in Lemma [2/6) 
Now we write this procedure as an algorithm to be used in the rest of the paper. 
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Algorithm 1 — Quasi GCD Algorithm 

Input: {{hi, • • •, hk }, ho, Y} where ho,h\,... ,hk £ IK n [y] and hi = hijY 3 for i = 

1) • • • > h, j 0 ,... ,d 1. 

Output: B = {7o, • • •, 7jvi}, where % = {{},{ h i,j, 1 < i < k,0 < j < d},{h 0 }}> T q = 
{W-Kb-- • > 9q,N 2 }>{5g,o}}(! < Q < Nx), such that 

Ah 

Zero(/ii,.. .,h k /ho) = (J Zero(^ 9 , 5 3il ,..., g q ,N 2 /g q ,o) 

5=0 

where 'L 0 = 0, deg('I' f/ . Y) > 0, and deg (g q ^, Y) = 0 for 1 < q < Nx, 0 < i < N 2 . 


Example 2.8 We use a simple example to explain the algorithm. Let the original polyno¬ 
mial system be {Y 2 + Y/Y}. First, we introduce a new variable Y\ and get an equivalence 
system {Y 2 + Y. Y\ Y — 1}. Second, we introduce a new variable Yq to make it homogeneous 
{Y 2 + YYq , Y\Y — Y^}. Finally, we introduce U,Uq,U\ and add YU + YqUq + Y\U\ to the 
homogeneous system. The matrix A corresponding to the homogeneous system is 


1 

0 

0 

0 

0 

0 

u 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

Uq 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Ux 

0 

0 

0 

0 

1 

0 

-1 

0 

0 

0 

u 

0 

Uq 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

u 

0 

Ux 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Uq 

0 

0 

Ux 

1 

1 

0 

0 

0 

0 

Uo 

0 

0 

u 

0 

0 

0 

0 

1 

1 

0 

0 

Ux 

0 

0 

0 

u 

0 

0 

0 

0 

0 

0 

-1 

0 

U\ 

0 

0 

0 

Uo 

0 

0 

1 

0 

1 

0 

0 

0 

0 

Ux 

Uq 

u 


A (nu m) 

is the submatrix of A formed by the first 6 columns, rank{A^ num <) = 6. According 
to Lemma \2. r /\ we must choose the first 6 columns and by calculating we find the submatrix 
formed by the first 9 columns and the last column is nonsingular, which is 


1 

0 

0 

0 

0 

0 

u 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

Uq 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Ux 

0 

0 

1 

0 

-1 

0 

0 

0 

u 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

u 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Uo 

Ux 

1 

1 

0 

0 

0 

0 

Uq 

0 

0 

0 

0 

0 

1 

1 

0 

0 

Ux 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

Ux 

0 

Uq 

0 

0 

1 

0 

1 

0 

0 

0 

0 

u 


det(A) = —Uf(U — Uq + U\). Substituting U\ =(),[/ = —1, Z7o = Y to U — U$ + U\, we 
obtain the polynomial —1 — Y and Zero(y 2 + Y/Y ) = Zero(E + 1) 

The components of Lemma 12 .6 1 mav be empty, as shown by the following example. 
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Example 2.9 Let h\ = xy + l,h -2 = x, and take y as the maximal variable. According to 
Lemma UTR it can be divided into two components Zero(l, x) and Zero(xy + 1, x/x). We can 
delete the first component. However, we cannot delete the second component Zero(xy+l,x/x) 
which is empty. The second component will be deleted in our main algorithm later when we 
continue our procedure to Zero(x/x). 

2.3 The decomposition algorithm 

We now give the main result about polynomial systems. 

Theorem 2.10 For a given polynomial system H = {h ±,..., hk} £ K n , deg(hj) < d for 
l < i < k, there is an algorithm to compute regular triangular sets A q = ..., ] 

which have the following properties: 

1. Zero (H) = U^ 1 Zero(sat(^4 g )) is an unmixed decomposition. 

2. The degrees of ’Lgp,..., \h Qy i q are less than d° n , N < k n d ncn+ ~, where c is a constant. 

n+2 

3. The running time of the algorithm can be bounded by a polynomial in k n and d nc 

Using the algorithm described below, we can calculate the regular triangular sets A q 
which satisfy the properties in Theorem 12.101 


Algorithm 2 — Algebraic Triangular Decomposition 
Input: {hi,...,hk}, where hi, ...,%£ K[xi,..., x n ]. 

Output: R, which is the set of (A q ,D q ) and A q = {ihgp,... Af q j q }. 'h q ,i,D q £ for 

l<g<A, l<i</g such that A q are regular triangular sets, ini(\l/ 9i j) / 0 on 1. 
any element of Zero (A q /D q ), and Zero(h \,..., hk) = U g Zero(sat(Ag)). 

Let T = {{}, {h u ..., h k }, {}}, § = {T}, R = {}. 

2. If § = 0 output R, else let T = {F, P, N} £ §, and § = S \ {T}. 

3. If |F| > k, go to step 2. 

4. If P = 0, add (F, OpeNlO t° R and go to step 2. 

5. Let x 7 = maxfcgp lv(/r), P = {h £ P 11 v(h) = x 7 }, P = P \ P. 

6. Let N = {/ £ N | lv(/) < x^}, H = U f- 

fen 

7. Apply Algorithm 1 to {P ,H,x 1 }, and let the output be ID. 

8. If D = 0, go to step 2, else let 71 = {W, U, V = }} £ D, B = D \ {71}. 

9. Let U = U U P, x v = maxj G u lv(/). 

10. If lv(u) < x v or U = 0, add {F U W, U, N U ¥} to S. 

11. If lv(u) > x q , write v = as a multivariate polynomial in x = ( x v ,... ,® 7 ) with 

coefficients in R[xi,... ,Xp_i]. Add {F U W,U,N U VU {/ Q }} to § for each a. Go to step 8. 


Example 2.11 A simple example is used to explain the algorithm. Let f = xyz + 1 ,g = 
x 2 + x, x < y < z. In step 5, x 7 = z and P = {/}. In step 7, applying Algorithm 1 to P, 
the output is Bi = {71} where 71 = {{xyz + 1}, {}, {xy}}. In step 9, we have U = {x 2 + x} ; 
x q = x. Since lv(xy) = y > x, we execute step 11 and add {{xyz + 1}, {x 2 + x}, {xy,x}} 
to § and go to step 2. Now we have T = {F,P, N}, where F = {xyz + 1}, P = {x 2 + x}, 
N = {xy, x}. In step 5, we have x 7 = x, P = {x 2 + x}. In step 6, we have N = {x}, 
H = x. Applying Algorithm 1 to {P, H, x}, the output is {{x + !},{},{}}• In step 10, we 








add {{xyz + 1, x + 1}, {}, {xy}} to §. In step 4, since P = 0, we add {{xyz + 1, x + 1}, {xy}} 
to M and output ®L Finally we have 

Zero (xyz + 1, x 2 + x) = Zero (xyz + 1, x + 1 /xy) = Zero (sat (xyz + 1, x + 1)). 

The purpose of step 11 is to add x to N. Otherwise, we will apply Algorithm 1 to {{x 2 + 
x},xy,x}, which does not satisfy the input condition of Algorithm 1 since x < y. 

Before proving Theorem 12.101 we first prove several lemmas. 

Lemma 2.12 Algorithm 2 terminates and each A q is a triangular set. 

Proof: By Lemma 12.61 after step 7, for any {W, U,V} G B, we have lv(p) < x 7 for any 
p G U. In other words, for any new {F,P, N} added to § in steps 10 and 11, the class of 
the polynomials in P will be decreased at least by one. Therefore, the algorithm terminates. 
Also, W is either empty or W = {p} and lv(p) = x 7 , which means A q is a triangular set for 
each q. □ 

Lemma 2.13 Omitting Step 3, Zero(/ii,... ,hk ) = 0^ = {Leio{A q /D q ), and ini(Tq i j) A 0 on 
any element of Zero(A q /D q ). 

Proof: To show Zero(hi,..., h k ) = (J Zero(A q /D q ), it suffices to show that the equality 

9 

Zero(hi, ...,h k ) = U {FiPiN}e §Zero(F U P/ JJ p) (8) 

peN 

always holds in the algorithm, and when P = 0 the algorithm returns the requires equa¬ 
tion. § is modified in steps 7, 10 and 11. In step 7, by Lemma 12.61 Zero(P /H) = 
U{w,u,{tj}}eiD>Zero({W UU/u). Clearly, after applying Algorithm 1, ([8]) remains valid when 
P and N are properly replaced as in steps 10 and 11. In step 1, a special substitution is 
performed. Let v = T,l a x. a . Then Zero(/u) is replaced by U Q Zero (/l a v). Since Zero(/u) = 
U a Zero (/l Q v), <[8j) is still valid after step 11. 

Now suppose (\hi,..., dq/M) is one component of the output. From the procedure of 
the algorithm, we know that this component is obtained in the following manner: 

Zero(/o,i,...,/ 0)fe (o)/M 0 ) —s- Zero^i, f lt i ,..., f hk w/Mi) 

->■ Zero(Ti, T 2 ,/ 2 ,i, ■ ■ ■,/ 2! fc( 2 )/MiM 2 ) (9) 

—)• ... 

-»■ Zero('Li,'L 2 ,...,VM 1 ...M t ) 

and M = Mi ... M t . Note that after applying Algorithm 1, Zero(Ti, /iq,..., f l ,fc ( i)/Ti) is a 
component of Zero(/o,i,..., / 0 fc (o))/Mo). If lv(u) < x v in step 10, Mi = T\. Otherwise, Mi 
is the multiplication of 7\ and a coefficient l a of Ti as shown in step 11. The component 
Zero(T 2 , / 2 ,i,..., / 2>fe ( 2 )/M 2 ) is obtained similarly from Zero(/ 2 j i,..., / 2jfc ( 2 )/S 2 ), where S 2 is 
the maximal factor of Mi satisfying lv(S , 2 ) < lv(/ 2j ) for all j. Continuing this procedure, 
we will obtain Q. It is obvious that 

Zero(/o,i,...,/ 0 )fe (o)/M 0 ) D Zero^i, /pi,..., / ljfe( i)/Mi) 

D Zero(Ti, T 2 ,/ 2 q,...,/ 2 jfc ( 2 )/MiM 2 ) (10) 

D ... 

D Zero(Ti,d' 2 ,...,d' t /Mi...M i ) 
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According to (1) of Lemma 12.61 we have ini('I'j) | Mi , so ini(\Eb) A 0 on any element of 
Zero(^i, □ 


Lemma 2.14 The triangular sets A q = ... ^t q ,i q } are regular and Zero(A q /D q ) = 

Zero(sat(^l (? )). 

Proof: Let Zero(4'i,..., 4 ’t/M) be a component of the output. According to the proof of 
Lemma 12.131 this component comes from procedure ®. Now we assume that lv(\Hi) = 
x kl 1 M\ G K fcl _i, 1 v(4' 2 ) = x k2 ,M 2 G ..., lv(4b) = x kt ,M t G K kt _i. 

According to Lemma [2.31 to show that A q is regular, it suffices to prove that inif+j) is not 
always zero on any irreducible component of sat(4'j + i,..., 'lb) for 1 < i < t — 1. We prove 
this by induction. First, supposing ini(\lL-i) is zero on an irreducible component of sat(4'j). 
Zero(d'f/Mt) is a component of Zero(/t_i i i,... , f t _ 1 k (t-i)/St~ 1 ) after applying Algorithm 1, 
where S t -1 is a factor of M t _ 1 . Obviously, Zero(4' i /AL 4 ) is not empty and Zero(4' t /M i ) = 
Zero (4b/ini (4b)) = Zero(sat(4b)) since \v(M t ) < lv(4/+ Since ini(^F t _i) is always zero 
on an irreducible component of sat(4b), there exists an r] kt = (£ 1 , ■ • • Ak t ) in Zero(4b/Mt) 
such that ini(^ t _i)(ry fct ) = 0. Since Zero(/i_i i i,..., f t _ lk{ t-i)/S t -i) D Zero(4b/M t ), r] kt G 
Zero(/ t _i.i,..., f t _ 1 k ( t -i)/St-i)■ If Zero(4b/M t ) is obtained from step 10, then r/ = rj kt 
is also in Zero(/ t _i ) i, ■ ■ ■, f t _i k (t-i)/M t -\). Otherwise, Zerois obtained from step 
11, 7j kt can be extended to a zero 7 / = rj kt _ 1 -i of Zero(/t_i i i,..., f t _ l k o-i)/M t -i), since 
St -1 is a coefficient of M t -1 • So in each case M^ t -i)(v)^ 0) but we have ini(4h_i)|M i _i, a 
contradiction. We have proved {\Ib_i, ^>t} is regular. We can prove in the same way that 
Mt-i is not always zero on any irreducible component of sat('Lj). According to Lemma 12.51 
we have Zero(4't_i, v h i /ini('h 4 )ini('& i _i)) = Zero(4't, ^t-i/Mt-iMt). The induction step can 
be proved similarly. □ 


Lemma 2.15 In Algorithm 2, the degree of the polynomials ..., are less than d° n 
and N < k n d nc , where c is a constant. The running time of the algorithm can be bounded 
by a polynomial in k n and d ncn . 

Proof: According to Lemma [2761 for given polynomials h\,...,h k G with deg(L,) < d, 
after applying Algorithm 1, we obtain no more than kd cn components, each component has 
no more than kd cn polynomials, the degrees of polynomials in these components are less 
than d c , and the running time of the algorithm can be bounded by a polynomial in k,d n . 
After applying Algorithm 1, the most complicated situation is that the maximal leading 
variable of the polynomials g q j is x n -\. Applying Algorithm 1 to these components, each 
component will be split to at most A:d m d c2 ^ Tl_1 ^ < kd° 3n components, each component has at 
most kd cn d c f n_1 ) < kd c n polynomials, and the degree of each polynomial is less than d c . 
This procedure will terminate in at most n steps. In step n, each component will be split 
to at most kd c " n components, each component has at most kd c " n polynomials, and each 
polynomial has degree less than d c ". Then in total, there are at most k n d c " +2n components, 
and the degree of the polynomials can be bounded by d° n . The running time of Algorithm 
2 can be bounded by a polynomial in k n , d cn+2n . □ 

Proof of Theorem \2. 1(A Omitting Step 3, the correctness of the theorem follows from 
Lemmas 12.121 12.131 12.141 and 12.151 It suffices to show that with step 3, the theorem is also 
correct. Suppose A k , k = 1,... , Nq are the extra regular triangular sets obtained by omitting 
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Step 3 and Ai,l = 1,... ,N are those obtained with Step 3. Then 

Zero(fai,..., h/.) = U^ 1 Zero(sat(.A/)) [^J U^ 1 Zero(sat(v4.fc)). 

From the condition |Fl > k in Step 3, we have |4fc| > k. By the dimension theorem proved 
in [J], d i rn (Z e r o (s a t (Afc))) < n — k. While by the affine dimension theorem [14^ p. 48], any 
component of Zero(/ii,..., h &) is of dimension no less than n — k. Thus, Zero(sat(4fc)) are 
redundant in the decomposition and can be deleted. □ 

3 Decomposition of ordinary differential polynomial systems 

In this section, a decomposition algorithm for ordinary differential polynomial systems will 
be given, which has an elementary worst case complexity bound. 

3.1 Basic definition and property 

Let IK be a field of characteristic zero in which an operation of differentiation is performable 
such that for any a, it £ I, 


(a + b)' = a' + b 1 , ( ab)' = ab' + ba!. 


Then we call IK a differential field. Let y\,... ,y n be differential indeterminates. We write 


the j-th derivative of y* as yf \ Let IK{yi,..., y n } = IK [y\ 3> ,i = 1,...,n; j G N] be the ring 
of differential polynomials in y\,, y n ■ 

Let / be a differential polynomial in IK{yi,... ,y n }- The class of / denoted by els(/), 

is the greatest p such that some y p ^ is present in /. If / G IK, then els(/) = 0. The order 

of / w.r.t yi, denoted by ord (f,yi), is the greatest j such that y\^ appears effectively in /. 

We write ord(/) = maxi<j< n ord(/, yi). If els(/) = i and ord(/, y*) = j then we call 

the leader of /, and we write it as ld(/) = yf\ We define y\^ = ld(/) > ld(y) = y^ if 

i > a or i = a,j > /3. We can write / as a univariant polynomial in its leader such that 

/ = ad(y ( / ] ) d + ... + do, and we call ad the initial of /, which is denoted by If. We call -^4r 

dy\ J 


.O') 


the separant of /, which is denoted by Sf. For f,g€ IK{yi,..., y n }, we say / is of higher 
rank than g, if one of the following conditions is satisfied 

1 . els(/) > els {g). 

2 . els(/) = els( 5 ) = p and ord(/,y p ) > ord (g,y p ). 

3. els (/) = cls(» = p, ord(/, y p ) = ord (g,y p ) = j, and deg (f,y ( p j) ) > deg (g,y { p ] ). 

Let cls(g) = p > 0. We say / is reduced w.r.t g, if ord(/, y p ) < ord(g, y p ) or ord(/, y p ) = 
ord (g,y p ) = j, and deg (f,y p j) ) < deg (g,y p j) ). 

For G K{yi ,... ,y n }, we use [fi,..., fk] to denote the differential ideal gener¬ 

ated by /],..., j )., which is the linear combination of /],..., ff, and their derivatives. 

A set T := {Ti,..., T r } of differential polynomials in IK{yi ,... ,y n } is called a triangular 
set, if cls(Tj) 7 ^ c\s(Tj) for i ^ j. Assuming that cls(Ti) < ■■■ < cls(T r ), we rename the 
variables as u\,... ,ut,y\, ■ ■ ■ ,y r such that r + t = n and ld(Tj) = ■ A differential 

polynomial / G IKjiti,..., Uf, yi,. ■ ■, y r } is said to be invertible w.r.t T if [/, Ti,..., T r \ FI 
IK{«i,..., ut} / {0}. T is called regular if It.- are invertible w.r.t to 77 —1 for 0 < i < r. T is 
called saturated if T is regular and S 7 ^ are invertible w.r.t to Tj for 1 < i < r. 
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Let T := {Ti,... , T r } be a triangular set. Denote Ij- = It x ■ • • It t and Sj- = ■ • • St t - 

Then the saturation ideal of T is 

dsat(T) = {/ € K{ yi ,..., y n } \3d £ N, s.t. (I T S T ) d f G [Tj,..., T r ]}. 

It is known that if T is saturated, then dsat('T) is an unmixed radical differential ideal [2|[3j. 

Lemma 3.1 Let T := {Ti,..., T r } be a triangular set in K{yi,..., y n }. Then T is saturated 
ifd-Ti and St, are not identically zero on all irreducible components o/dsat(7i_i) and dsat(7I), 
for 1 < i < r, respectively. 

Proof: This lemma can be proved similar to Lemma 12.31 □ 

Lemma 3.2 Let T := {Ti,..., T r } be a saturated triangular set in lK{t/i,..., y n }. If M € 
K{yi,... ,y n } is not identically zero on all irreducible components o/dsat(T), then we have 
Zero (T/Lt'S't) = Zero(T /M Ip St) = Zero(dsat(T)). 

Proof: This lemma can be proved similar to Lemma 12.41 □ 


3.2 A squarefree quai GCD algorithm 

In order to decompose differential polynomial systems, we need to modify Lemma 12.61 In 

Lemma [2Tl for given polynomials ho, h\,... ,hk € K[xi,..., x n , Y], deg(hj) < d, we can write 
d -1 

hi(i > 0) as hi = hijY J , and divide the whole space as K n = (J Uij [J{x G K n | h t j (x) = 
j= 0 ij 

0, V 1 < i < k and 0 < j < d — 1}, where 


U^j — Zero(hi j d_i,... , hi$, ■ ■ •, ^2,o 5 • • 


for 1 < % < k, 0 < j < d — 1. We write hij = 

0</3<j 

polynomial system becomes 


■ j hi t d— i, ■ ■ • , hi t j-i-i/hij) 

hi,pY^. Then on Uij, the original 


— A - l — • • • — hjs — 0, /iq ^ 0. 

We add a step here to divide (flTT) into the following polynomial systems: 


( 11 ) 


hij = h 
hij = h 

hij — h 


i+l — • 

• • — hk 

o 

O 

II 

i+l = ■ ■ 

II 

'c§ ^ 
ii 

i+l = ■ ■ 

II 

*ci ^ 

H 


d -h±^ o 

dY T 

d 2 hi j 

■ =0 ’h 0 —^yl 0 


dY 2 

d j ~ l h 




OY'J 


-i 


0, h 0 / 0. 


( 12 ) 


Since = hij, and hij ^ 0 on Uij , we actually have 0. Then the zero set of (TTTT) 

equals to the union of the zero sets of (1121) . Now we continues to introduce new variables as 
in Lemma [2.6l to make the polynomial systems homogenous. After this modification, Lemma 
12.61 becomes the following form. 
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Lemma 3.3 Given polynomials ho, hi,..., hk G K n [T], deg(hi) < d, and hi = J2j =o 

for 0 < i < k, we may compute g q j, G K n , Tg G K n [Y] \ K n for 1 < q < Ni,0 < t < N 2 such 

that: 


Ni 

Zero(/ii,.. .,h k /ho) = Zero(^ g ,p gi i,.. ■,3 9 ,iV 2 /Sg,o)UZero({/i i j, 1 < i < k,0 < j < d}/h 0 ) 
g =1 

which has the following properties: 

1. We have ini('Fg) | g qfi , and S,i. q / 0 on any element of Zero(d> q , g q j, g qi N 2 /g q ,o)- 
2- degx lv ..,x n ,y(^ , <2)? d e gxi,...,x n (5?,t) A V(d)-,Ni,N 2 < kV(d n ). 

3. The running time of the algorithm can be bounded by a polynomial in k,d n . 

Proof: For property 1, we need only to prove that S^ q / 0 on any element of Zero(\kg, g q \,..., 
g q .N 2 )- Since we divide (fill) into the union of f[T2]) . each component of the output, for example 
('hi, <71 ,... ,5x2/50 )) comes from one of (fT2l) . Without loss of generality, suppose it is the 
first one in (1121) . Then we have 


~ on' ' 

Zero (hij, h l+{ ,..., h k /h 0 -^r) D Zero(Ti, ST ...,g N2 /g 0 ). 

If Sq fl vanishes on (£1 ,... ,f n Al) G Zero(\hi, </i,... ,gN 2 /go), then rj must be a multiple root 
of Ti when substituting (x\,... ,x n ) by (£1 ,..., f n ). According to Lemma [2.71 rj is also 
a multiple root of the homogeneous equation system of (1121) after introduce new variables 

Yi,Yq, which means ...,£ n , rf) = 0, a contradiction. Property 1 has been proved. 

Property 2 comes from LemmaE© We now prove property 3. According to the procedure 
of this algorithm, the origin system has been divided into no more than kd subsystems Hij 
in (©. For each H, j, we divide it into no more than d subsystems in (1121) . and each system 
has no more than k + d polynomials, and the degree of these polynomials are bounded by 
2d. The related matrix A has , 2 rows, where 

D = ( (7; - !)) + To < 6 d. 

l</<min{ 2 ,fc—2+1} 


The degree of the elements in A are bounded by 2d, so the degree of P s in Q and A s in 
© are bounded by 2dC]- )+2 < d(6d + l)(6d + 2). Since P s are linearly independent, we have 

s < (d(6d+ l)(6d + 2)) n ). For each W s in ©, we divided it into Wp in © and l is no more 
than the degree of A. So in total, we have kd 2 (d(6d + l)(6d + 2))^” +1 ^) components and 
Ni < kV(d n ). According to the above proof, it is obvious that the degree of each polynomial 
in these components is no more than the degree of P s and A, so is bounded by V(d). The 
polynomials g q j. come from three parts. The first part is the polynomials in U t j and whose 
number is bounded by kd ; the second part is the coefficient of P s when taken as polynomials 
in U, Uq, U\ and so the number is bounded by (d(6d + l)(6d + 2)) 2n ); the third part is the 
coefficients of Ep and so the number is bounded by (d(6d + 1)(6 d + 2)) n+1 ). Therefore, 
N 2 < kV(d n ). □ 

Now we write this theorem as an algorithm. We only give the input and output of this 
algorithm, since the procedure of this algorithm has been described above. 
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Algorithm 3 — Squarefree Quasi GCD 

Input: {{hi,...,h k },{h 0 },{xi,...,x n },Y} where h 0 ,h 1 ,...,h k € K[xi ,... ,x n ,Y], 

deg (hi) < d. 

Output: D = ,T Nl }, where To = {{},{hij,l < i < k,0 < j < d - l},{h 0 }}, 

Tq = {{’Lg}, {g q , 1 , • • • ,g q ,N 2 }, {3g,o}}(l < q < Ni), which satisfy the conditions 
in Lemma 13.31 


3.3 The algorithm 

We now give the main result for differential polynomial systems. 

Theorem 3.4 Let hi,...,hk € K{yi,... ,y n }, where deg (hi) < d and ord (hi) < R for 
l < i < k. There is an algorithm to compute saturated triangular sets A q := ..., ^ q ,i q 

which have the following properties: 

1. Zero (It i,..., hjf) = U^L 1 Zero(sat(Ag)). 

2. We have deg(® gii ) < d c2 " H ; ord(®,,i) < 2 n R, and N < k 2nR df' RRn . 

3. The running time of this algorithm can be bounded by a polynomial in k 2nR d c2 RRn . 

We will give an algorithm to produce those saturated triangular sets in the theorem. 
Before giving the main algorithm, two sub-algorithms will be given. The first one is the 
partial remainder mm- 


Algorithm 4 — DPM Algorithm 


Input: {{go}, {/i, •••,/*},{/()}}, where g 0 , f 0 , ■ ■ ■, fk € K{yi ,..., y n }, ord (fi,y a ) < r 

for 0 < i < k, and Id (go) = Va t > 1. 

Output: {{g 0 , /l, ■ ■ ■, fk}, {fo *Sg 0 1} where ord (A, y a ) _< r - t for 0 < i < k such that 
Zero (g 0 , /i, ..., fk/foS go ) = Zero (g 0 , /i , - - -, fk/foS go ). 


1. For i = 0,..., k, 

1.1. A = 

1.2. If ord(/j,g Q ) < r — t, goto step 1. 

1.3. Let ord(A, y a ) = r* and g^ i ~ r+t ' > = S' 90 y£ r<) - H n . 

1.4. Replace y£ r, '' ) in A by and multiply by (S go ) des ^ i,y ° ,x) \ and let A* be the new 
differential polynomial. Goto step 1.2. 

2. Output {{g 0 , /i,..., Afc}, {Ao<Sgo}} • 


Lemma 3.5 07] / Use the notations in Algorithm 4 and assume deg (fj) < d, deg(go) < d, 
ord (/i,y 7 ) < R for 0 < i < k, 1 < 7 < n. Then, we have the following bounds: ord(/j,g 7 ) < 
R + t, deg(fj) < V(d, t) for any 0 < j < k, 1 < 7 < n. 

Next, we describe a splitting subroutine from nn. Let g € K{y 1 ,... ,y n }. For a € 
{1, •••,«}, let ord(g, y a ) = r, g = daiVaVa^ ■ ■ ■ Va ] ) a , a = (a 0 , a u ..., a r ), (y a y^ ... 
Ua ) ) a =Va° ... {ya ) ) ar . Denote coeff(g, y a ) to be set of g ija . For G C K{yi, ..., y n }, denote 

coeff (G,y a ) = U geG coeff(g, y a ). 

We have the following split algorithm. 
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Algorithm 5 — SPLIT Algorithm 

Input: {G,y a }, where G = {g u ...,gi}c K{y u ... ,y n }. 

Output: B = {%,...,Tn} where % = (coeff(G, y a ), 0), % = ({A,i,..., A/J, {{nAy}) 

dy a 

such that ord(hj.i, y a ) = > 0 and 

dh 

Zero(<7i, ...,g t ) = U^ 1 Zero(/i ii i,... ,h ih / —^A-) U Zero(co eS(G,y a )). (13) 

dy% l) 

1. Let § = {{gi,... ,g t }}, B = 0. 

2. If § = 0, return B; else let F G § and S = § \ {F}. 

3. If V/ G F, ord(/, y a ) = 0, then add (F, 0) to B. Go to step 2. 

4. Let / G F such that ord(/, y a ) = t > 0, deg(/, y$) = d. Set F = F \ {/}. 

5. Let / = J2j=o l d(ya ) ) d and A = d ^{)y A = 1, • • • , d. 

6. Let B = B(J {({/} U F, {A}), ({A,/} U F, {/ 2 }),..., ({/*£_!,. • •, /} U F, {A})}, and S = 
S(J{F U {Zq, ii, • • •, Id}}- Go to step 2. 


Note that the order and degree of the difference polynomials in the output are smaller 
than or equal to that of g j in the input. We now give the decomposition algorithm. 


Algorithm 6 — Differential Triangular Decomposition 
Input: {hi, ...,h k }e K{yi ,... ,y n }. 

Output: { A q , Dq}, 1 = 1,... , N, where A q = ... 4 fqj } satisfies the conditions in 

Theorem m 

1. Let T = {{}, {hi,... h k }, {}}, § = {T}, R = {}. 

2. If § = 0, output R, else let T = {F,F, N} G § and § = § \ {T}. 

3. If P = 0, add (F. O p gnP) to R and go to step 2. 

4. Let yi 7) = max feeP Id (h), P = {h € P | ld(h) = yi 7) }, P = P \ P. 

5. Let N = {/ G N | ld(/) < y^}, H = n /eS /• 

6. Apply Algorithm 3 to {P, H, vars(P U {H}) \ {yi 7 ^}, ya' 1 }, the output is B. 

7. If B = 0, go to step 2, else for A = {W = {V’},U,¥ = {u}} G B, B = B\{A}, U = UUP. 

8. If U = 0, add {FUW,U,NU ¥} to R and go to step 7. 

9. Apply Algorithm 5 to {U, y a }, the output is Bi. 

10. If Bi = 0, go to step 7, else let C = (T, 0) G Bi and Bi = Bi \ {C}. 

11. If 0 A 0, assume © = {-%r}. Applying Algorithm 4 to {{<?}, W U (r \ {#}),¥}, the 

_ &ya 

output is {W,¥}. Add {F, W, N U¥} to §. Goto step 10. 

12. If 0 = 0, let y(^ = max/j g r ld(/i), y^ = ld(u). 

13. If y A' ) < y^ add {FUW, r,¥U N} to §. Goto step 10. 

14. Let x = {Ay | deg {y,y^') > 0 and y^ e > ye} and write v as a multivariate polynomial 
in x: v = S/^x 0 . Add {F U W, T, N U ¥ U {le}} to § for each 9. Go to step 10. 


We use two examples to illustrate Algorithm 6. 


15 












Example 3.6 Note that Algorithm 6 can also be used to algebraic polynomial systems and 
return a radical decomposition. Let f = xy 2 with x < y. Using Algorithm 2 to f , we 
obtain two components {x} and {y 2 }, where the second one is not radical. In step 6 of 
Algorithm 6, when applying Algorithm 3 to {{/}, {}, y}, the system {/ = 0} is first split into 
{xy 2 = 0, 2 xy 0), {xy 2 = 2 xy = 0, x 0), and {x = 0} and then returns 0, {{?/}, {}, {x}} ; 
and {{},{s},{}}. Finally, we obtain the decomposition Zero(/) = Zero(x) UZero(y). 

Example 3.7 Let f = y ' 2 — xy 2 , x < y. In step 4, we have ya ' 1 = y', P = {/}. In Step 
5, N = 0 and H = 1. In Step 6 , Algorithm 3 is applied to {P ,H,y'}. P is first split into 
two components {y ' 2 — xy 2 = 0, 2 y' 0} and { y ' 2 — xy 2 = 2 y' = 0}. The output of the first 
component is {{y ' 2 — 4 xy}, {}, {xy 2 }} and the output of the second one is {{y'}, {xy 2 }, {}}. 
In Step 8 , Co = {{y ' 2 — xy 2 }, {}, {xy 2 }} will be put into § and eventually be added to R. 
In Step 9, we will handle {{y'}, {xy 2 }, {}}. Applying Algorithm 5 to U = {xy 2 } and 
Ua = V, the output Bi consists of C\ = ({xy 2 }, {2xy}), C 2 = ({xy 2 , 2xy}, {2x}), and C 3 = 

«M,{})- 

C\ is handled in Step 11. Algorithm 4 is applied to {{xy 2 }, {y'}, {2xy}} and returns 
{{xy 2 , x'y 2 }, {2xy}}. Finally, C 4 = {{}, {xy 2 , x'y 2 }, {2xy}} is added to S. 

Co is handled in Step 11. Algorithm 4 is applied to {{2xy}, {y', xy 2 }, {2x}} and returns 
{{2xy,2x'y, xy 2 }, {2x}}. Finally, C§ = {{}, {2xy, xy 2 ,2xy, x'y}, {2x}} is added to §. 

C 3 is handled in Steps 12 and 13. Co = {{y'}, {2x}, {}} is added to S. 

For 64 , in Step 6 , Algorithm 3 is applied to {{xy 2 , x'y 2 }, 1, {2xy}} and returns the empty 
set. We omit the computing procedures for Co and Cq. The algorithm give the decomposition 
Zero(/) = Zero(sat(/)) U Zero(y',x) U Zero(y). 

Now we prove Theorem 13.41 with the following lemmas. 

Lemma 3.8 Algorithm 6 terminates, Zero(/ri,..., hjf) = UgZero(^4 g /D g ), andlq, qi , Sq, qi 7 ^ 
0 on any element of Zero(A q /D q ). 

Proof: The algorithm has three loops, starting at Step 2 , Step 7, and Step 10, respectively. 
We need only to show that the loop starting at Step 2 will terminate. Let {Fi,Pi,Ni} be a 
component added to § in this loop and yf = maXpgPj ld(p). Then, we have either y & c < y2 
which means that the algorithm terminates. 

Zero(/ii,..., h]f) = {J Zero(4' (?] i,... y V q j q /D q ) can be proved similar to Lemma 12.131 In 
g 

the proof, we also need the equalities in Lemma 13.31 Lemma 13.51 and (1131 . 

We now show that A q is a triangular set. It suffices to show that for any {F,P, N} £ S, 
maxpgp cls(p) < maxqgp cls(y). New polynomials are added to P in Steps 8 , 13, and 14. In 
Step 8 , since U = 0, this is indeed the case. In Steps 13 and 14, we have 0 = 0 which 
means that y a and its derivatives do not appear in T. Hence max pG r cls(y) < a and A q is a 
triangular set for any q. 

Finally, if (Ti,..., Tt/M) is one component of the output, then according to the algo¬ 
rithm it comes from a procedure like ([9]) and (HOI) . In the algebraic case, from one step to 
the next step in Q, Algorithm 1 is used one time. In the differential case, from one step to 
the next step in ([9]), Algorithm 3 is used many times. For instance, the procedure to obtain 
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'I'i is as follows: 


Zero(/ 0 ,i,... ,/ 0)fc (o)/Mo) ->• Zero(^i, hip, ■ ■ ■, h 1)t (i)/Mi) 

->• Zero(3i )0 ,...,5i;(i)/S'iMi) 

->• Zero(V> 2 , ^ 2,1 > ■ ■ ■, h um /M 2 S 1 M 1 ) 

-> ... (14) 

-t Zero(5 Si0 , • • • ,g Sjl (a)/M s+1 • • • Si Mi) 

-1 Zero^g+ij/ig+i,!,..., /i s+1)t(s+ i)/S s+ i¥ s+ i • • • Si Mi) 

where 'Pi = V’s+i and {/i s+ i,i, • • •, h s+lt(s+ i)} = ■ ■ ■, f ljk m} in ©■ (V’l, hip, • ■ ■, 

h lt (i)/M\) is a component of (/op,..., f 0 k { 0 )/M$) after using Algorithm 3, so J^jMi and 
S- l p 1 A 0 on any element of Zero(^i, hip,..., h 1 fc (i)/Mi) by Lemma HUfl (gig ...., g 1 gi)/S\M\) 
is a component obtained from (ip i, hip,..., h 1 k (i)/M{) by Algorithms 5 and 6 in Steps 9 and 
11. So we have Zero(<7po, • • •, pp/SiMi) C Zero^i, hip,..., h l fe (i)/Mi). The procedure 
is repeated until cls(ip s+ i) > cls(g s _|_ip) for all j and Ti is obtained. Then, we have 

Zero(/ O p,..., f o k (o)/Mo) D Zero(^i, hip,..., h htW /Mi) 

D Zeio(g lfi ,...,g ll(1) /S 1 M 1 ) 

2 ... (15) 

D Zero(g S)0 ,.. .,g s y a) /S s • • ■ SiMi) 

D Zero(^ s+ i, h s+1 p,..., h s+ljt (»+i) /M 8+1 • • • Si Mi) 

By Lemma [Mil Iq l |M s+ i and Sy 1 A 0 on any element of Zero(^ s+ i, h s+ ip,..., h s+l t (s+i) /M.s+i 
■ ■ ■ Si Mi). The lemma is proved. □ 

Lemma 3.9 In Algorithm 6, A q := IP^p, • • • T qj are saturated triangular sets and 
Zero(dsat(^lg)) = Zero(^ g p,... 'P q ,i q /D q ). 

Proof: Using Lemmas 13.11 and 13.21 instead of Lemmas 12.31 and 12.41 the proof of this lemma is 
the same with that of Lemma 12.141 □ 

The following lemma gives the complexity part of Theorem 13.41 

2 n 

Lemma 3.10 In Algorithm 6, the degree of is less than d c , the order of is 

less than 2 n R, N < k 2nR d c2 RRn . The running time of this algorithm can be bounded by a 

polynomial in k 2nR d c2 RRn . 

Proof: For given differential polynomials hi,..., h k € K{yi,..., y n }, since ord(h,, yj ) < R for 
1 < i < k, 1 < j < n, there are at most Rn variables when applying Algorithm 3. Consider 
the most complicated case where yf ; is the maximal leader. After applying Lemma 
[3731 to hi,..., hfc, there are at most kd cRn components, each component has at most kd cRn 
differential polynomials, and the degree of each polynomial is no more than d c . After applying 
Algorithm 5, each component will be split into at most d cR components, and each component 
has at most kd cRn+cR polynomials. After applying Lemma 13.5( the number of the components 
and the number of polynomials in each component do not change, and the degrees of the 
polynomials are less than d 2c . The most complicated case occurs when the order of y n 

/ d _o\ 

decreases by one and the maximal leader becomes yg . Continuing this procedure until 
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the main variable becomes y n -\. There are no more than k R d cRRn components in total, each 

component has no more than kd c Rn polynomials, the degree of these polynomials are less 

than d c , and the order of these polynomials are less than 2 R. Repeating this procedure 

to 2 / 1 ,... ,y n - 1 , finally we obtain no more than k 2 ’ lR d c2 RRn components, the degree of the 

2 n 

polynomials are bounded by d c , and the order of these polynomials are less than 2 n R. □ 


4 Summary 

Two triangular decomposition algorithms are given in this paper. For a set of polynomials 
F = {/i,..., f s } in K[xi,..., x n ], we can compute regular triangular sets 71,..., T r such that 
Zero (F) = Uj Zero (sat (%)) which gives an unmixed decomposition for the solution set of F = 
0. We also show that the complexity of the algorithm is double exponential in n. For a set 
of ordinary differential polynomials F = {/i,..., f s } in K{yi,..., y n }, we can give a similar 
decomposition Zero(F) = UjZero(dsat(7i)), where dsat(7j) are radical differential ideals and 
the complexity is triple exponential. This seems to be the first triangular decomposition 
algorithm for differential polynomial systems with elementary computation complexity. 
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